Is the distant globular cluster Pal 14 in a deep-freeze? 
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ABSTRACT 

We investigate the velocity dispersion of Pal 14, an outer Milky- Way globular cluster at 
Galactocentric distance of 71 kpc with a very low stellar density (central density 0.1-0.2 Mq/pc^). 
Due to this low stellar density the binary population of Pal 14 is likely to be close to the primordial 
binary population. 

Artificial clusters are generated wit h the observed pro perties of Pal 14 and the velocity dispersion 
within these clusters is measured as Uordi et al.l ( 20091) have done with 17 observed stars of Pal 14. 
We discuss the effect of the binary population on thes e measurements an d find that the small 
velocity dispersion of 0.38 km/s which has been found bv lJordi et al.l (|2009l ) would imply a binary 
fraction of less than 0.1, even though from the stellar density of Pal 14 we would expect a binary 
fraction of more than 0.5. We also discuss the effect of mass segregation on the velocity dispersion 
as possible explanation for this discrepancy, but find that it would increase the velocity dispersion 
further. 

Thus, either Pal 14 has a very unusual stellar population and its birth process was significantly 
different than we see in today's star forming regions, or the binary population is regular and we 
would have to correct the observed 0.38 km/s for binarity. In this case the true velocity dispersion 
of Pal 14 would be much smaller than this value and the cluster would have to be considered as 
"kinematically frigid", thereby possibly posing a challenge for Newtonian dynamics but in the 
opposite sense to MOND. 
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Stellar populations with densities less than 
about 10^ stars/pc^ have a binary fraction. 
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as in the Taurus-Auriga groups (density about 
10 stars/ pc'^) where the fraction is near 100 per 
cent (e.g. iKroupa fc Bouvien 120031 ). But even in 
high-density star-forming regions, such as the 1 
Myr old Orion Nebula cluster which has a den- 
sity of 10''' stars/pc'^ ([Hillenbrand &: Hartmann 
[1998 ), the binary fraction is near 50 per cent 



of at least 50 per cent. For example, in the so- 
lar neighbourhood (density 1 s tar/pc'^) /; 



(jKohler et al.l 120061 ). Stellar-dynamical models 



0.6 (JDuquennov fc Mavon Il991[ ) or in OB as- 



sociations (density < 0.1 star/pc"^) fbin > 0.7 



([Kouwenhoven et all |2007| ). In star- forming re 



gions the binary fractions can be even higher, such 
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have demonstrated that the Orion Nebula cluster 
must be expanding and that consequently further 
destruction of binary systeins has mostly ceased 
(|Kroupa. Aarseth. fc Hurlevll200ll ). Only in glob- 
ular clusters, which are dense (> 10^ stars/pc"^) 
and which have a long dynamical history over a 
Hubble time, are binary f ractions observed to be 
less than 10-20 per cent (JHut et al.l [19921) . But 
this is probably understandable as a result of the 
destruction of binary systems, first in the very 



dense initial config uration which these objects typ- 
ically formed with (^Marks. Kro upa. fc Baumgardt 



20081 : iMarks fc Kr oupa 20 la), followed by long- 



term binary-star burning. 

Pal 14 is an old globular cluster with a mass of 
about 10^ Mq, a (three-dimensional) half-light ra- 
dius of 34 pc, and is located at adistaiiceof 
71 kpc from the Galactic centre ( Hilkeii l2006t 
IJordi et al.ll2009l) . It has a low density (0.1-0.2 
stars/pc'^) and so conformity with known stel- 
lar populations would imply it to have fbin ~ 0.5. 
Measurements of fbin would thus allow testing the 
dependence of the binary population on the phys- 
ical conditions during star formation, as Pal 14 is 
a low-metallicity population that f ormed nearly a 
Hubble tini e ago (IJordi et al.ll2009l) . Furthermore, 
Ijordi et al.l (J2009|) suggest that Pal 14 may be sig- 
nificantly mass segregated. How mass segregation 
can affect a velocity dispersion measurement has 
not been investigated, though. 

However, Pal 14 is also interesting for testing 
gravitational theories. That binaries can affect 
the dynamical mass estimate of stellar populations 
has already been studied in the context of dwarf- 
sphero idal satellite galaxies by iHargreaves et al 



(11996) and has recently beeii re-a ddressed by 



Kouwenhoven fc de GriisI (|2008l ) and iGieles et al 



(|2010l ) for star clusters in the Milky Way disk. 
The measured (one-dimensional) velocity disper- 
sion of Pal 14 is (0.64 ±0.15) km/s or (0.38 ±0.12) 
km/s, depending on whether the outlying "star 
15" of a sample investigated bv ljordi et al.l (J2009I ) 
is included or not, respectively. Now, if Pal 14 has 
a binary fraction of 50 per cent or more, then the 
velocity dispersion measured from a single spec- 
troscopic snapshot may be significantly contam- 
inated, i.e. inflated, because of the binary-star 
orbital motions of the tracer stars. A potential 
problem arising here is that the dynamical velocity 
dispersion of the cluster must be about 0.5 km/s 
for Newtonian dynamics to hold. So there is not 
much room for a binary-star contribution to the 
measured velocity dispersion. Alternatively, if the 
binary population is normal (i.e. canonical), then 
the true (binary-corrected) velocity dispersion of 
the virialised cluster may be much smaller than 
the Newtonian value. If this were the case, then a 
Newtonian crisis would emerge, such that Pal 14 



would be "kinematically frigid" , a situation which 
is not expected from current theoretical dynamics. 
A third possibility is given by the unreliability of 
low-number statistics in such systems, i.e. a sam- 
ple of about 17 stars may simply be insufficient to 
determine the velocity dispersion of such a stellar 
system. 

Thus, either Pal 14 is in virial equilibrium, such 
that its stellar population would be non-canonical 
(a binary fraction much lower than all other known 
populations). Or, the cluster has a canonical stel- 
lar population, but is then kinematically frigid, 
in violation of Newtonian dynamics. This would 
not be remedied by MOND, because according 
to MOND the cluster ought to have an even 
larger velocity dispersion th an the Newtonian 
value (JBaumgardt et al.l 120051 ). The third expla- 
nation, of low-number statistics, would imply that 
velocity dispersion measurements of low-velocity 
stellar systems would have to be taken with much 
more caution than is done in practice nowadays. 

With this contribution we assess the expected 
velocity dispersion of Pal 14 assuming various 
binary fractions and Newtonian dynamics to be 
valid (Section 2). These numerical experiments 
are compared to the observed velocities in the 
cluster in Section 3. The conclusions are con- 
tained in Section 4. 

2. Velocity dispersions 

For the purpose of testing the observed velocity 
dispersion of Pal 14, we create artificial star clus- 
ters with the observed properties of Pal 14 and 
measure the velocity dispersion in these clusters. 
Therefore we use the code McLusteeQ (Kiipper 
et al., in prep.), an open-source project developed 
at AlfA Bonn to set up clusters for A^-body com- 
putations, especially for NbodyG (Aarscth 20031). 
We use this software to conveniently generate 300 
renditions of Pal 14 with various properties for 
taking out mock observat ions and testing the find- 
ings of Ijordi et al.l ( 20091 ) on their conclusiveness. 



The parameters f or the artifi cial clu s ters w e take 
from IJordi et al.l (|2009[ ) and iHilkeiJ (|2006[ ). The 
most important quantity in this respect is the 
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total mass of the cluster as it should be closely 
linked to the velocity dispersion in the assumed 
case of the cluster being in virial equilibrium. Us- 
ing star counts, Ijordi et alj (J2009I ) find a mass for 
Pal 14 between 6000 - 12000 Mq without taking 
into account compact remnants, and depending 
on the assumed mass function within the clus- 
ter. Assuming addit ional 15-20 per cent ma ss in 
stellar remnants fflabringhausen et al.ll2009l and 
Dabringhausen private communication) yields a 
mass range for Pal 14 of 7000 - 14000 Mq. We 
therefore concentrate on the two cases of 7000 Mq 
and 14000 Mq, respectively. 

As the density distribution we choose a Plummer 
sphere with a half-mass radius of 32 pc, which 
i s the observ ed half-light radius of Pal 14 of 1.'28 
(lHilkeiil2006h at the as sumed distance of (71 ±1.3) 
kpc (jjordi et al.l 120091 ). Setting the half-light ra- 
dius equal to the half-m ass radius a s sumes that 
mass follows hght, just as I Jordi et al.l (|2009l) have 
done in their investigation, i.e. that the cluster 
is not significantly mass segregated. Moreover, 
the infinite Plummer distribution is cut off at the 
Jacobi radius (- 128/156 pc for 7000/14000 Mq) 
assuming a Galactic circular velocity of 220 km/s 
at the Galactocentric radius of Pal 14. 

As the mass function of the cluster stars we use the 



canonical IMF (JKroupal 120011 ) . but cut it off at a 
maximum stellar mass of 1.0 Mq as all stars above 
this mass should have died by stellar evolution at 



the e xpected age of Pal 14 of 11.5 Gyr (|Jordi et al 



20091 ). Even though ljordi et al.l (|2009l ) find a shal 
lower slope in the range of 0.5 Mq to 0.8 Mq, we 
argue that the actual slope of the mass function 
is not very important in this respect as the veloc- 
ities of the cluster stars get drawn independently 
of each other and we, in the end, only observe 
stars of mas s abo ve 0.7 Mq to be consistent with 
(|2009[ ). 



Jordi et al 



For investigating the effect of a realistic binary 

population on the velocity disper sion, we set up 

the bi nary population following the lDuquennov fc Mavor 
( 19911 ) period distribution for field stars and a 
thermal eccentricity distribution. The binary 
orbital planes are distributed randomly, as are 
the orbital phases. We reject binaries with 
a semi-major axis below 100 Rq, though, as 



these binaries may have been destroyed in a 
common-envelope p hase of the binary components 
(|Zorotovid lin prep.l ). Moreover, we use random 
pairing for the binary components since we are 
dealing with an evolved population of low-mass 
stars. We vary the binary fraction, fbim from to 
1 and generate ten clusters for each binary frac- 
tion to gain better statistics. 

In addition, we set up clusters with the above 
properties but being mass segr egated. Therefore , 
we use the procedure defined in lSubr et al.l (|2008[ ) 
whic h is implemented t hrough Plumix in McLus- 
TER. ISubr et al.l (J2008l ) defines the degree of mass 



segregation through a single parameter, S, which 
can vary from (not segregated) to 1 (completely 
mass segregated) . To keep the number of models 
to a minimum, we concentrate on the cluster with 
7000 Mq without binaries, and see how mass seg- 
regation affects its velocity dispersion. Again we 
generate 10 clusters for each value of S. 

The line-of-sight velocity dispersion, a, we mea- 
sure from our clusters by calculating 



^2 - v^. 



(2) 



where v is the mean velocity in the sample, and 
v'^ is the mean squared velocity. We do this sep- 
arately in three different directions and take each 
as an independent measurement. 



In Fig. [T] we show the mean velocity dispersions 
for all stars above 0.7 Mq for ten different binary 
fractions and for both cluster masses. The error 
bars give the standard deviation of the different 
cluster renditions from the mean. The clusters 
without binaries have a very little spread about 
the mean whereas already a binary fraction of 0.1 
introduces such a large scatter that the two mass 
groups of clusters overlap within their error bars. 
At a binary fraction of 0.5 the two mass groups are 
indistinguishable as their velocity dispersion gets 
completely dominated by the binary population. 
At a binary fraction of 1.0 the velocity dispersion 
is about four times larger than in the case without 
binaries. 

In Fig. [2] we show the velocity dispersion when we 
draw 17 stars o ut of the tot a l sam ple of stars above 
0.7 Mq, just as I Jordi et all (|2009l ) have done. The 



data points show the mean of 3 000 000 renditions 
out of the ten clusters for each binary fraction. 
Again we took the hne-of-sight velocity dispersions 
independently along three different directions for 
each rendition. Since the final distribution of ve- 
locity dispersions for each binary fraction does 
not follow a Gaussian distribution but is rather 
asymmetric, the error bars show 68 per cent of all 
renditions below the mean and 68 per cent of all 
renditions above the mean. 

We rejected stars, though, when their velocities 
were off the mean by more than 2.5 km/s, i.e. only 
took stars within a velocity window of Av = 5 
km/s into accoun t. This is similar to the way 
Jordi et al.l (|2009l ) get their velocity dispersion of 
Pal 14, and is the usual practise to discriminate 
between cluster members and non-members. We 
see that the mean measured velocity dispersion 
can increase by about 50 per cent compared to 
the velocity dispersion of all stars above 0.7 M© if 
we only have 17 stars in our sample. This effect is 
small for low binary fractions but gets significant 
above a binary fraction of about 0.3. Moreover, 
the scatter in our velocity dispersion measurement 
is huge. At fbin = 0.5 we can measure values be- 
tween 1 km/s and 3 km/s, depending on which 
stars we take into our sample. At a binary frac- 
tion of 1.0 the mean measured velocity dispersion 
is about 7 times larger than the Newtonian veloc- 
ity dispersion without binaries. 

In Fig. [3] we show the same experiment but with 
a velocity window of Av — 10 km/s. The ef- 
fect described above gets more significant. The 
mean value of the measured velocity dispersion in- 
creases and also the scatter grows. At fbin = 0.5 
we can get values between 1 km/s and 4 km/s. 
At a binary fraction of 1.0 the measurements of 
a small sub-sample and the measurements of all 
stars barely agree within the error bars any more. 
Note that the mean value can be as large as 10 
times the true Newtonian velocity dispersion with- 
out binaries. Moreover, from comparing Fig. [2] to 
Fig.|3]we see that the measured velocity dispersion 
depends strongly on the velocity window which we 
allow. We tend to increasingly overestimate the 
true velocity dispersion by increasing the size of 
the allowed velocity window. 



In Fig. |4] the velocity dispersion measurements 
are shown for the mass segregated clusters. The 
plot shows the velocity dispersion measurements 
for a cluster of 7000 Mq without binaries, as seen 
in Fig. [T][3] but for a mass segregation degree, S, 
varying from to 0.95. From the figure we see 
that mass segregation increases the velocity dis- 
persion further. The segregated cluster (5 = 0.95) 
has a 20 per cent higher velocity dispersion than 
the unsegregated case {S — 0). This is due to the 
fact that through mass segregation the heaviest 
stars (which we observe for the velocity disper- 
sion) move to the cluster centre, and that stars on 
average move at higher velocities when they are 
in the cluster centre. When a sample of observed 
stars was concentrated on the centre of the cluster 
and was not well distributed over the cluster, we 
would get the same effect. 

3. Discussion 



Jordi et al.l (|2009l ) have deter mined the velocity 



disper sion of Pal 14 according to iPrvor fc Mevlan 
(J1993I ). i.e. have made a maximum- likelihood esti- 
mation, which is, of course, necessary as unlike in 
our samples the measurement errors of their radial 
velocities are all different. From their sample of 17 
stars they find one star (star 15) to lie 3 cr off the 
mean value, and thus split their investigation into 
two parts: one with taking star 15 into account 
and one without taking star 15 into account. For 
the sample with star 15 they find a velocity dis- 
persion of (0.64±0.15) km/s, and without star 15 
they find (0.38±0.12) km/s. 



Based on a Kolmogorov-Smirnov test, iJordi et al. 
( 2009[ ) argue that star 15 is most likely not a reg- 
ular member of Pal 14, i.e. a foreground contam- 
ination or part of a long-period binary system. 
Thus, at the bottom line they favour the lower 
value of 0.38 km/s. 

From Fig. [T] we see that such a low velocity dis- 
persion would imply that the binary fraction of 
Pal 14 was less than 0.1. With taking star 15 into 
account the velocity dispersion of Pal 14 would 
imply a binary fraction of less than 0.2. 

In Fig. [2] the velocity dispersion derived from a 
sub-sample of 17 stars within a velocity interval of 
Av = 5 km/s shows that the lower value of 0.38 



km/s would be consistent within the error bars 
with values less than fbin = 0.2, while the higher 
value of 0.64 km/s would be consistent with values 
less than fbin = 0.4. Allowing for Av — 10 km/s 
reduces the latter to fbin less than 0.3 (Fig. [2|). 

From Fig. |4] we see that mass segregation does 
not help in understanding the low observational 
value of Pal 14's velocity dispersion, since mass 
segregation tends to increse the observed velocity 
dispersion of a cluster by up to 20 per cent com- 
pared to the unsegregated case. On the contrary, 
since Pal 14 is su pposed to be mass segregated 
(jjordi et al.l 120091 ) . its velocity dispersion may be 
even inflated. Our previous estimates on Pal 14's 
binary fraction therefore have to be taken as up- 
per limits. 



4. Conclusions 

From our test of the observed velocity disper- 
sion of the Milky- Way globular cluster Pal 14, we 
have seen that the binary fraction within Pal 14 
has to be less than 0.2 in the case of star 15 not be- 
ing considered a member, in order to be consistent 
with a velocity dispersion as low as (0.38±0.12) 
km/s. With taking star 15 into account the max- 
imum binary fraction consistent with the obser- 
vational uncertainties is about 0.4. This poses a 
number of questions on the nature of Pal 14. 

As a first explanation of these findings we may 
assume that the binary fraction within Pal 14 is 
indeed as low as found. But, as stated above, the 
density within Pal 14 is as low as 0.1-0.2 stars/pc"^. 
The effect of disruption of binaries due to dynam- 
ical stellar evolution within its age of about 11.5 
Gyr is therefore negligible and we even may see 
here the p rimordial binary population of a globu- 
lar cluster ( Hasani Zonoozi et al.l2010l ). Thus, the 
formation of Pal 14 must have been significantly 
different from what is observed in all other star 
forming sites today. Or Pal 14 must have under- 
gone a very dense and violent phase in which most 
of the binaries were burned, but as we can see in 
the ONC today, this is very unlikely to happen for 
a sufficiently long time span that the binary frac- 
tion drops below 0.5. Furthermore, such a scenario 
would give rise to the question how Pal 14 could 



have expanded that much, as with a half-light ra- 
dius of about 34 pc it is one of the most expanded 
globular clusters of the Milky Way today. Recent 
numerical studies moreover show that this expan- 
sion is very unlikely to be of pure dynamical origin, 
since expansion takes place on a relaxation time 
scale and the relaxa tion time of Pal 14 is of or - 
der of a Hubble time ( Hasani Zonoozi et al.ll2010l l. 



A second option would be that the binary popu- 
lation is normal, i.e. above 0.5, and that thus the 
observed velocity dispersion has to be corrected 
for the effect of binaries. As we have seen, this 
can be as much as a factor of 10 in the case of 17 
stars drawn from a cluster with a binary fraction 
of 1.0. Since the lower mass li mit of Pal 14 is de- 
termined to be about 7000 M© (jJordi et al.ll2009l ) 
this would imply that the true velocity dispersion 
of Pal 14 is much lower than the Newtonian pre- 
diction, i.e. Pal 14 is "kinematically frigid" . This 
would be inconsistent with our understanding of 
Newtonian gravity and could neither be explained 
by considering MOND to be valid in Pal 14. 

Moreover, we found that mass segregation in- 
creases the observed velocity dispersion of a clus- 
ter even further, and thus cannot explain the low 
velocity dispersion of Pal 14. The observed, flat- 
tened mass function of Pal 14 on the other hand 
sugges ts that Pal 14 is significantly mass segre- 
gated (jJordi et al.ll2009l ). Thus, its unsegregated, 
i.e. for the effect of mass segregation corrected, 
velocity dispersi on may be even low er than the val- 
ues reported by Ijordi et al.l (|2009h , which would 
worsen the problem. 

Thus, the present state of knowledge on Pal 14 
is that either its binary fraction is highly abnor- 
mal, given its low density, or that it is significantly 
sub-virial. The former possibility would imply a 
non-standard star-formation event which formed 
Pal 14, while the latter indicates a problem un- 
derstanding the dynamics of Pal 14. 

In any case, this investigation has shown that 
one has to be cautious with low velocity disper- 
sions derived from small samples of cluster stars. 
The observed value from such a sub-sample tends 
to be significantly larger than the true velocity 
dispersion. This effect gets larger with increas- 



ing size of the allowed velocity window about the 
mean radial velocity. This could be especially im- 
portant in the ou ter parts of star clusters (e.g. 
Scarpa et al.ll2Q07 ). as was also recently shown by 
Kijpper et alJ(|2010[) . 



More line-of-sight velocity measurements, e.g. as 
a 2nd epoch spectroscopic snapshot, may help 
to reduce the statistical uncertainties and im- 
prove the significance of findings on whether or 
not Newtonian dynami cs is valid in Pal 14 (see 
also iGentile et al]l2010l ) . 
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Fig. 1. — Mean velocity dispersions of the cluster 
renditions with 7000 M© and 14000 Mq as deter- 
mined from all stars above 0.7 M0 for ten differ- 
ent binary fractions, fun- The error bars give the 
standard deviation from the mean. The velocity 
dispersion can increase by a factor of 5 for large 
binary fractions, compared to the binary free case. 
Above fbin = 0.5 the binaries dominate the veloc- 
ity dispersion such that the two cluster groups get 
indistinguis hable. Also shown are the measured 
results from Ijordi et al.l (|2009l ) with and without 
taking into account star 15 (see text). Their veloc- 
ity dispersion measurements allow a binary frac- 
tion of less than 0.1, or 0.2 taking into account 
star 15, respectively. 
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Fig. 2. — Same as Fig. [T] but additionally the 
velocity dispersions are shown which are derived 
from sub-sets of 17 stars above 0.7 M©. For this 
measurement all stars lying outside a velocity in- 
terval about the mean value of Aw — 5 km/s 
were rejected. The data points show the mean of 
3 000 000 renditions, and the error bars show the 
values in whic h 68 per cent of a ll renditions lie. 
The results of I Jordi et al.l (|2009() allow a binary 
fraction of less than about 0.1, or 0.3 taking into 
account star 15, respectively. 
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Fig. 3. — Same as Fig. [2] but for a vefocity in- 
terval of Av — 10 km/s. The effects from Fig. [5] 
get even more significant. The mean derived ve- 
focity dispersion from our artificial clusters can be 
as large as 10 times the Newtonian velocity dis- 
persion in the case without binaries, and is always 
significantly larger than the true velocity disper- 
sion within the cluster. Even though the spread 
about this mean value is large, at fbin = 1.0 the 
velocity dispersion of the sub-sample and of all 
stars barely agree within the error bars. 
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Fig. 4. — Mean velocity dispersions of the cluster 
renditions with 7000 M© as determined from all 
stars above 0.7 M© (as in Fig.[T]), as well as derived 
from sub-sets of 17 stars above 0.7 M© in a veloc- 
ity intervall of Av = 5km/s (Fig. ^ and Av = 10 
km/s (Fig. 131), respectively, but here for ten dif- 
fer ent degrees of ma ss segregation, S, as defined 
by ISubr et all (|2008() . For increasing mass segre- 
gation the velocity dispersion rises as heavy stars 
are preferentially located in the cluster centre and 
hence move at a higher velocity on average, thus 
mass segregation cannot e xplain the low vel ocity 
dispersion measurement of Ijordi et al.l ( 20091 ). 



